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Introduction 

« This report is divided into two major sections. The 

first section considers in some detail the problem of 
atmospheric interference on remotely sensed data. A program 
is described which allows the effects of atmospheric 
absorption and emission by carbon dioxide, water vapor, 
methane, nitrous oxide, carbon monoxide, and ozone to be 
determined so that correction techniques can be devised. 

The program covers the spectral region from 0.69 to 20um 
and divides the atmosphere into up to thirty layers of equal 
optical mass. The theory needed to understand the program, 
the program description and operating instructions, and 
comparisons to test data are given in the first section of this 
report. A complete description of the program including a 
prbgram listing, input card formats and sample data runs, 
and samples of the data used in the program are given in ' 
Appendix A of this report. 

The second section of this report gives an analysis of 
ocean truth data and subsequent laboratory analyses which 
were performed in support of this study. Unfortunately, a 
coordinated analysis of Sky lab acquired data could not be. 
performed with associated ocean truth data as was planned 
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for our two different test sites. For test site 1 which 
was located in the Eastern Atlantic Ocean adjacent to the 
northwest coast of Africa the ocean truth data for this 
investigation was to be acquired and supplied to us by the 
Principal Investigator of EPN 518. Unfortunately, due to 
the shifting track of the Skylab spacecraft, EPN 518 was 
unable to acquire ocean truth in conjunction with a 
Skylab overpass and, thus, EPN 518 was unable to supply our 
investigation with the required ocean truth data. The 
second test site for this investigation was located in the 
Gulf of Mexico in the near coastal waters off the coast of 
Freeport, Texas. The actual location of this site shifted 
during the course of Skylab from the coastal waters of Freeport, 
Texas to those of Matagorda Bay, Texas, and finally to the 
offshore waters off the coast of Matagorda Bay, Texas. 

A total of four ocean truth expeditions were conducted 
beneath the Skylab spacetrack on 5/31/73, 6/5/73, 8/8/73, and 
finally on 9/16/73. Of these four expeditions, Skylab 
acquired data on only two occasions, namely 6/5/73 and 8/8/73. 

Of these two data sets, the data on 6/5/73 was badly 
contaminated by solar glitter'*', and the data from 8/8/73- 
was contaminated by approximately 95% cloud cover. Thus, 
since no Skylab data could be acquired with the required 
concomitant ocean truth data, the results of section two of 
this report are limited to a presentation and analysis of the 
ocean truth data obtained. 

The solar elevation angle on 6/5/73 was 82® and the 
maximum solar elevation angle required by the Investigation 
ReQuirem ent Document was 65® . 


1 


T 


ATMOSPHERE TRANSMISSION COMPUTER PROGRAM CP^ 

By D. E. Pitts, T. L. Barnett, C. L. Kerb, 

Walter Hanby,* and Alyce E. Dillinger* 

Lyndon B. Johnson Space Center 

SUMMARY 

All Earth resoxirces remote-sensing techniques are affected, to some 
degree., by the atmosphere lying between the sensor and the test site. The com- 
puter program described herein allows the effects of atmospheric absorption and 
emission by carbon dioxide, water vapor, methane, nitrous oxide, carbon monox- 
ide, and ozone to be assessed so that correction techniques can be devised and 
tables for atmospheric correction algorithms can be calculated. The absence of 
carbon monoxide data on the data tape, however, limits calculations to the 
remaining five gases. The program incorporates data from wavelengths, of O .69 
to approximately 20 micrometers, generated by R. Calfee of the National Oceanic 
and Atmospheric Administration, but can incorporate other line-by-line data in 
the proper format. The atmosphere is divided into layers (530), each of equal 
mass j and is defined by significant levels of either radiosonde data or a 
model atmosphere. Both upwelling and downwelling emissions are calculated to 
enable simulation of both terrestrial and spaceborne measurements. The program 
is designed so that punched cards can be produced for use in other computer 
programs. Sample runs for both radiosonde data and laboratory data are pro- 
vided in the appendixes. ' 

INTRODUCTION 


The Objectives of the NASA Earth Observations Program are to determine the 
performance capabilities of various sensors, to identify signature criteria of 
resources, to develop new sensors and systems, to devise analysis systems and 
procedures, and to demonstrate applications using these tools. The laboratory 
usually affords the best testing environment for a sensor, but the type of tar- 
get, the conditions of the signal path, and other testing parameters are lim- 
ited. Because the' laboratory environment is usually restrictive, a successful 
laboratory test of the sensor is necessary but not sufficient to ensure proper 
.operation of the sensor in the real environment. Therefore, much of the test- 
ing is performed in 'an environment that duplicates the conditions under which 
the instrument is expected to operate. The success of the testing under such 
conditions requires that the data concerning the environment between the 
instrument platform and the test site be as accurate as possible. 


^This report was produced under Skylab EREP contracts EPN-502, 58A, & 586. 
•Lockheed Electronics Company, Inc. 
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The magnitude of the effect on the signal caused by interaction between 
the atmosphere and the signal depends on the type of sensor used, the wave- 
length employed, and the meteorological conditions prevailing at the time of 
the measurement. Furthermore, the atmosphere can cause large deleterious 
effects in some analysis systems (e.g., maximum-likelihood pattern recognition 
computer programs); thus, the analysis system as well as the sensor and elec- 
tronics must be evaluated. Therefore, the output from the CP computer program 
has been arranged so that pianched cards can be produced for use in other com- 
puter programs to predict the effect of the atmosphere or to correct the data 
for atmospheric effects for a data collection and analysis system. 

The computer program described herein allows the calculation of atmos- 
pheric transmission due to carbon dioxide, water vapor, methane, nitrous oxide, 
ozone, and carbon monoxide. Because no data exist on the data tape for the 
last of these, however, carbon monoxide calculations are precluded. The pro- 
gram incorporates data from wavelengths of O .69 to approximately 20 micrometers 
but can incorporate other line-by-line data in the proper format. 

Appreciation is expressed to Dr. R. F. Calfee, who assisted the authors in 
using his model and data, and to David Anding, who kindly ran his band model 
for the Salem, Illinois, test case. 


SYMBOLS 


A 

BOUND 




full width at half maximum for triangular slit function, cm ^ 

distance from a line at which the wings of the line must be con- 
sidered (usually approximately 12 cm~^), cmT^ 

constants determined from the boundary conditions 


CAYBUR 
CONN = 


dummy variable 


RO 


10 / 

c speed of light, 3.0 x 10 cm/sec 

DELV the increment in cm”^ for calculations of upwelling and downwelling 

A 

intensity, - — ; also distance the center of triangular slit is 
moved to make the next intensity calculation. 


DV 


integration step within the triangular slit, cm 


-1 




E„(x) 


E" 

g 




Ov 




Kourganoff function, 


-1 


lower rotational energy level of the transition, cm 
acceleration due to gravity 

2 

black-body monochromatic intensity at T and v, W/(cm •sr*vjm) 

2 

monochromatic intensity of Sun or other source, W/ (cm 'sr'ym) 

monochromatic upward-welling radiation (intensity or radiance), 
2 

W/(cm •sr*ym) 


monochromatic downward-welling radiation (intensity or radiance). 


W/(cm •sr'vmi) 


degraded intensity, W/(cm •sr'ym), averaged over frequency 


K 


atmospheric layer index 

dtunray variable representing atmospheric layers 
monochromatic mass absorption coefficient. 


(molecules /cm ) 


-1 


V - In 2 

“■0 ■ Tt 


k 

£ 

H 

M* 

m 


Boltzmann constant 


number of integration steps from v. - A to v + A, 

2A 11 

which is equal "to ~ 


mass of molecule, g 

molecvilar weight of the atmosphere, g/(g*mole) 
one interval in frequency 


3 


NN 

n 

P 


P. 

1 


number of equal- mass layers in the atmosphere 
monochromatic refractive index of the atmosphere 
atmospheric pressure, atm 

ambient atmospheric pressure for layer i, atm 


!• 

ii' 

^ I . 
n . 


0 

R 

RO 


S 

So 

SAY 

T 


1 atmosphere pressure, atm (lOl 325 K/m^, 1013.25 mbar) 
radius of planet 

universal gas constant based on the carbon-12 atomic weight scale 
in ergs/(K g*mole) 

variable factor in equation (19) » exponent which is a function of 
the molecular species 

—1 —2 
line intensity at P and T, cm /(molecule»cm ) 

—1 "2 

line intensity at and P^; cm /(g*cm ) for water and 

cm”^/ (molecule* cm) for carbon dioxide and other gases 

dummy variable 

temperature of the atmosphere as a function of z, K 
base temperature for spectral lines = 296 K for Calfee data 


temperature of the target, K 


transmission for ith layer 


i ,mono 


U 


transmission for all line-by-line calculations 


dummy variable of optical depth 


•/. 


optical mass, | p(z)dz, molecules/cm^ 


U. 


water- vapor optical mass, molecules/cm^ 


Wl. 

1 


water- vapor partial pressure, atm 


u 


En(x) 


Kourganoff function. 


1 

fo 


dy 


E" 

e 


I, (T) 
bv 


Ov 






lower rotational energy level of the transition, cm ^ 
acceleration due to gravity 

2 

black-body 'monochromatic intensity at T and v, W/(cm *sr'ijm) 

2 

monochromatic intensity of Sun. or other source, W/(cm *sr»ym) 

monochromatic upward-welling radiation (intensity or radiance), 

2 

W/(cm •sr'ym) 

monochromatic downward-welling radiation (intensity or radiance), 
2 

W/(cm •sr'um) 

2 

degraded intensity, W/(cm *sr*ym), averaged over frequency 


i 

K 


atmospheric layer index 

dximmy variable representing atmospheric layers 
monochromatic mass absorption coefficient. 


(molecules /cm ) 


-1 


ir L- ST 2 
0 ■ a. > IT 


k 

I 

M 

M* 

m 


Boltzmann constant ' 

number of integration steps from v - A to v + A, 

2A ^ ^ 

which is equal to ^ 

mass of molecule, g 

molecular weight of the atmosphere, g/(g ‘mole) 
one interval in frequency 


i 


W. 


weighting function for downwelling intensity (unitless), derivative 
of transmission with respect to altitude times Az. 


W 


weighting function for upwelling intensity (unitless), derivative of 
transmission with respect to altitude times Az. 


w 


X = 


variable factor in equation (l8) 
V - In 2 


X = -(t - t) 

. -- 

Y = -—-Vin 2 


altitude, cm 


a 


Lorentz half width at P and T, cm 


-1 




2kT In 2 


Az. 

1 


Me 


Lorentz half width at Tq and P^, cm 


height of atmospheric layer i 


-1 


Ap^ zenith extent of Sun or other source 

A<^^ azimutiial extent of Sun or other source 


gray-body emissivity assumed between and 


monochromatic emissivity of target 


6 

V 


zenith angle, rad 
cos 6 (unitless) 

position of Sun or other source (unitless) 


5 


i 




cos 0^ , vhere 0^ = angle through ith layer 


frequency, cm 


-1 


frequency of line center, cm 


-1 


frequency at beginning of calculation interval, cm 


-1 


P 

% 

o(v) 

a»(v) 

^Ov . 


h 

a 


frequency at end of calculation interval, cm”^ 

3 

constituent gas density as a function of z, molecules/cm 
monochromatic reflectivity of target 

normalized instrument function 
triangular instrument fvinction 

monochromatic optical depth of the entire atmosphere (dimensionless) 


monochromatic optical depth 


■ f 


p dz (dimensionless) 


azimuthal angle , rad 

azimuthal angle of Sun or other source, rad 
plateau function 


6 
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THEORY 


In many remote-sensing applications through the Earth atmosphere in 
which scattering is not an important phenomenon and local thermodynamic equi- 
lihrium can be approximated because of the lack of rapid chemical reactions, 
the equations for radiative heat transfer can be written in the manner of Love 
(ref. l), as indicated by the following equations and diagram. 


^4^ = -<(v“) ^ 






2 Top of the atmosphere 




where p = cos 0 


0 =: zenith angle 


= monochromatic upward-welling radiation (intensity or radiance) 


= monochromatic downward-welling radiation (intensity or radiance) 

.z 

= optical depth, K^p dz 


= optical depth of the entire atmosphere 


* monochromatic mass absorption coefficient 


7 


II 

p = constituent gas density as a function of altitude z 
^ = azimuthal angle 


n *= refractive index 



= black-hody intensity at atmospheric temperature 


T 


V = frequency 

The solutions to equations (l) and (2) for the case in which atmospheric 
temperature is a function of altitude require the use of an integrating factor. 




-T 


C^e + n 


t-T 
.X V 

V 


y 




(3) 


:(%-) = 




y bv 


I, . (t)dt 


ik) 


where and are constants determined from the boundary conditions and 

t is a dummy variable of optical depth. Assuming that the target is a diffuse 

reflector p having an emissivity e and a temperature T. and located at 
V VI 

t^ = 0, that Pg ” ^2^^2^ ~ ^ atmosphere 

(t^ *= and that the Sun or other source is illuminating the top of the 

atmosphere with radiance incident over the area represented by Ay and 

A()i at y^ and (j)^, respectively, the boundary conditions are as follows. 


, V p f2Trfl 


(5) 




(6) 
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THEORY 


In many remote-sensing applications through the Earth atmosphere in 
which scattering is not an important phenomenon and local thermodynamic equi- 
lihrium can he approximated because of the lack of rapid chemical reactions, 
the equations for radiative heat transfer can be written in the manner of Love 
(■ref. l), as indicated by the following equations and diagram. 



( 1 ) 




( 2 ) 


2 Top of the atmosphere ^ _ 


Ov 




Vr r rtTf T mT iT r r rrrrTrr r T T rr j t ,. = 0 


where y = cos 0 

6 = zenith angle 


= monochromatic upward-welling radiation (intensity or radiance) 

= monochromatic downward-welling radiation (intensity or radiance) 


Ty = optical depth, ^ K^p dz 


* optical depth of the entire atmosphere 


» monochromatic mass absorption coefficient 


7 


/3 


p s= constituent gas density as a function of altitude z 
= azimuthal angle 
n ~ refractive index 

I^^(T) = black-hody intensity at atmospheric temperature T 
V = frequency 

The solutions to equations (l) and (2) for the case in vhich atmospheric 
temperature is a function of altitude require the use of an integrating factor. 



where and are constants determined from the boundary conditions and 

t is a dummy variable of optical depth. Assuming that the target is a diffuse 
reflector having an emissivity and a temperature and located at 

= 0, that pg = 0 and ^ atmosphere 

(t^ *= and that the Sun or other source is illuminating the top of the 

atmosphere with radiance incident over the area represented by Ap and 

A(j>- at and respectively, the boundary conditions are as follows. 



(0,p) 


. V p r2-n ri 


(5) 
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'r 

where, when the plateau function = 0, (J> > ^ ^ 2^* 

4> < <|>1 - and y < “ (i.e., the function is a rectangular box in the 

middle of the (4 »m) plane); f2 = 1 otherwise. 

Solving for and in equations (3) and (4) by using equations (5) 

and (6) gives the following solutions. 




The two most important cases of remote sensing through an atmosphere are 
covered by equations (T) and (8): sensing upwelling radiance (eq. (T)) 

reflected and/or emitted from a target using downward-looking sensors and 
sensing downWelling radiance (eq. (8)) from targets. In both cases, the equa- 
tions are general and account for both monochromatic absorption and reradia- 
■tion by the atmosphere for targets at any altitude in or above the atmosphere 
(arbitrarily called t = 0) and for sensors in or above the atmosphere. 


l^v^bv(^l) TI 


2ti 


..1 


-T 


Ov 


Jiloe 


f'ov ^ ■ 

V' 

j— f 

y dy’ d<{)'>e ^ + n^ 

3 

/ 

•^0 


e ^ 


■I^^^(t)dt (7) 






T-t 

r ov — 

1, (t)dt 

U bv 


(8) 
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Equation (7) can be simplified somewhat by integrating over the azimuthal 
direction. 


T 11 


-T 




~ Ip A(^e ^ 
u O^v ^ 


-T 


Ov 


e ^ y» dy’ 


, Ay • 


-T.. /•! /•T 


+ 2p 


-T /•-!• 



2 y 
n e 

V 

Jr\ J 


Ov T-t 

y' 


e (t)dy' dt 


bv 


J /.T , ;^T 

r V* 

— 

0 ^ 


■\^(t)dt 


• (9) 


By taking the limit at small Ay for the upper radiation source and introduc- 
ing the exponential integral or Kourganoff function (ref. 2) 


En(x) = f y” 


-X 


-2 y , 
e dy 


(10) 
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where x - -(t - t), equation (9) then can be written 


-T 


-T 


V -T 




V 


Ov 




-r , >‘1 


-T yT 


+ 2p n^e ^ 
V 


Ov 


Eg(t - T)l^^(t)dt 


r\^ 

2 

Jn 


(t )dt 


(11) 


■, I 
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In the case of most natural targets, the value of is approximately 

0.9 or higher for wavelengths of 1+ to 20 micrometers and thus causes to 

be small. Thus, for remote- sensing measurements, the reflectance of atmos- 
pheric downwelling flux usually is negligible except when sensing is being 
attempted in or very near strongly absorbing bands from low-altitude platforms 
or when sensing is done at a very large angle from nadir, for which p is 
much larger. Thus, the third term in equation (ll) 


2p 


r'^ov 

AM = 

•'o 


E2(t - t)l^^(t)dt 


can usually be ignored. The second term in equation (ll) 


-T ”^0V 
V — 




11 


n 


is usually small as well, but is larger than the third term since the source 
considered is often the Sun. The second term shoxild always be a strong candi' 
date for inclusion at wavelengths <4 micrometers and shoiold be included for 
specular reflectances at any wavelength. 

For most purposes, it is therefore possible to simplify equation (ll) to 


-T 




= e 


L (T\e + 
V bvV 1/ 


2 

4 . S- 
h 


-T t-T 
^ V V 


e I, (t)dt 
bv 


( 12 ) 
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COMPUTER PROGRAM 


Program Development 


The terms in equations (8) and (12) are calculated in finite difference 
form by assiiming a refractive index of unity. The optical depth t in finite 
difference form is ^ 


NN 


r = / K p. Az. 


(13) 


i=l 


where NN is the mmiber of equal-mass layers i (usually NN t 10 ). The 
monochromatic intensity transmitted through the atmosphere (where j is a 

dumny variable representing- atmospheric layers and c is gray-body emissivity 
assumed between and Vg is 


. NN 

j=i 


tVi^) 


( 14 ) 
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i 


The upward-emitted intensity, designated in the printout (appendixes A and B) 
as "INT UP TOTAL” (for the NN-layer printout) is 



and the downward-emitted intensity, designated in the printout (appendixes A 
and B) as ”INT DOWN TOTAL” (for the NN-layer printout), is 




(16) 


Thus, the terms in equations (8) and (12) can he calculated for the Lorentz 
broadening of the line shape, where 



i/here k is the Boltzmann constant, T^ = 296 K, and = 1 atmosphere = 

1013.25 10^ N?m^ (1013.25 millibars); r is 0.62, 0,58, 0.5, 0.5, 0.5, and 

0.5, and w is 1.5, 1*0, 1.5, 1.5, 1.0, and 1.5, for water, carbon dioxide, 
ozone, nitrous oxide, carbon monoxide, and methane, respectively; E” is the 
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lower rotational energy level of the transition; is the line intensity in 

-1-2 u 

cm /(molecule* cm ); is the Lorentz half width; and is the freq.uency 

at a line center. These last four parameters are obtained from a magnetic tape 
containing data for carbon dioxide, water, nitrous oxide, ozone, and methane 
(ref. 3)« These data, a sample of which is shown in appendix C, are obtainable 
from the authors of this report on request. Other data can be used just as 
easily with the program, provided that the format and parametric units are the 
same or that statement 151 of the main program car. be changed to accommodate 
other data tape formats (ref. 4). The relationships between the units for 
optical mass U, for Kq(v), and for S indicated in the list of symbols in 

this document and the older units are given in table I, taken from reference 5* 

• Because the absorption at any frequency v results not only from lines 
near that frequency but also from the wings of nearby lines, the absorption 
coefficient must include contributions from these sources. This procedure is 
accomplished by summing all values of for ± BOUND, where BOUND is 

the distance from a line at which the wings of the line must be considered. 




*=Vq+BOUND 


’-Vq-BOUND V 


(v - V') 


2 u. 2 

+ a , 
v’ 


( 20 ) 


For low pressures, the predominant broadening effect is caused by the 
motion of the molecules (Doppler broadening). To account for this effect and 
for the transition between Lorentz and Doppler broadening, the Voigt profile 
(refs. 6 and 7) is used to calculate the absorption coefficient (K^) when 

|- < 0.25, X < 25.0, and Y < 5.0. 

0 



J 


] 



5^^ 



c = speed of light 


M s= mass of the molecule 

For very low pressures, equation (2l) approaches Doppler broadening; for higher 
pressures, it approaches Lorentz broadening. 

Because remote sensing is accomplished with instruments of finite band- 
width, equations (lU) to (l6) are integrated over a triangular instrument 
function a*(v) = A - jv^ - vj, where - A i v + A, as indicated by the 

following diagram, in which A ♦ is one-half the width of the slit triangle 
base, DV is the integration step with the triangular slit, and DELV is the 
step between printouts. 



For example, to find the degraded intensity I _ through the atmosphere 

(ref. 3) 



rVi+A 

I I o*(v 

Jv.-A ^ 

1 

rv.+A 

I ^ a*(v)d 
Jv,-A 


•V .+A 


a*(v)dv 


( 22 ) 
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where the normalized instrument function d(v) is defined by 


o(v) = 




/. 


v.+A 

^ o*(v)dv 


g*(v) 

.2 


( 23 ) 


v^-A 


and this step is performed for each Vq value from to Vg so that 


-Vj+A 

= I I^o(v 

•'v.-A 

1 


)dv 


(2U) 


Thus, the average intensity transmitted upward through the atmosphere plus that 
emitted for one position of the filter function at will be 




Ave 


m=£ NN 

m=l j=l 


-/k p, AzJ 

V "j ^ V 


‘J 


Avn 


- m=jl NN NN 

, ¥■ 

m*l i=l J=i 


"j 


( 25 ) 
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For a target of intensity Iq in or above the atmosphere, the average inten- 
sity transmitted through the atmosphere plus that emitted downward for one 
position of the filter function at will be 



where, NN = number of layers chosen for the atmosphere, i = one layer’ in the 

2A 

atmosphere, i = number of integration steps from - A to + A, 
and m = one interval in frequency. 

Thus, equations (25) and (26) describe the degraded intensities and 

• ^0 

I , which are upwelling and downwelling intensities, respectively, for a gray 

0 

body with emissivity e and temperature T^. 




(27) 



T -T 
V 




e 


y 


+ 



T-t 



by 


(t)dt 



( 28 ) 


The two terms in each equation are calculated in finite difference form 
and are printed and punched out by the computer program. To enable the 
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calculation of equations (25) and (26), two computer programs were integrated. 
One of these (ref. 3) is used to calculate a monochromatic atmospheric trans- 
mission spectrum for one layer of the atmosphere and then to degrade this spec- 
trum to the desired wavelength resolution as it might he observed with a 
spectrometer. The other program (ref. 6) affords a self-consistent method for 
calculating the state of the atmosphere on a refracted path through a spherical 
atmosphere, between two specified points, on the basis of radiosonde data or 
model atmosphere data (ref. 9). 


The combination of the two programs, using a layered atmosphere 

2 

(530 layers), produces high-resolution transmission and emission over any path 
by using realistic model atmospheres that can then be degraded to the appro- 
priate resolution. The model can be used to correct upward-, sideward-, or 
downward-looking instruments of either high or low spectral resolution in the 

frequency range of 0.7^2 to 1.UU15 >: 10^ cm The data include 15 250 lines 
of carbon dioxide, water vapor, ozone, nitrous oxide, and methane. Moreover, 
the model includes all significant levels of water vapor and as many as 
30 layers of differing concentrations of carbon dioxide, ozone, nitrous oxide, 
carbon monoxide, and methane. Furthermore, this computer program is not 
restricted to the Earth atmosphere. Any atmosphere for which any of the five 
constituent gases is appropriate may be studied by changing the molecular 
weight of the gas M*, the radius of the planet R, the acceleration of gravity 

6, and CONN, which equals where RO is the universal gas constant, to 

the appropriate planetary values in the main program data statement. 


Another calculation added was that of weighting functions for downwelling 
intensity (which is the derivative of transmission with respect to altitude 
times Az^) 




X-i. 

u 

4=1 



(29) 


2 " -1 
Resolution should not be pushed beyond 2 or 3 cm because Calfee com' 

presses the lines within each 1-cm interval into one or two equivalent 
"compressed" lines for each species. 

i 
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and for upwelling intensity 





J=i-1 



(30) 


These functions are required for inputs to inversion routines that derive 
temperature structure from the upwelling radiance in the iH- to 15-micrometer 
wavelength region and water-vapor specific humidity structure from the upwell- 
ing radiance at wavelengths of 6.3 to 7*0 micrometers and a knowledge of the 
temperature structure. Chahine (ref. 10), Barnett (ref. 11), and Smith 
(ref. 12) have shown that it is necessary to have weighting functions for a 
class of atmospheric conditions near the solution to ensure rapid convergence. 


Continuvun 

The wave number region between 800 and 1250 per centimeter encompasses the 
continuum absorption due to water vapor. Calfee has used the continuum absorp- 
tion coefficient given by Burch (ref. 13) for this region. First, he uses the 

2 

optical mass of water vapor (U^ in molec\iles/cm ) to calculate the partial 
pressure due to water vapor (Wl^) for the ith layer. 


Wl. = ~ 


UfTPi 


^ Az. * 7.2h9 X 10 
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(31) 


PH^ = -Wl^ + 0.005 * (P. + Wl^^ 


(32) 


CAYBUR = -0.0665335 x lo"^^ + 0.3721108 x l0*^^v^ 


-0.5971389 X 10"^^ V + 0.31131+85 X 10“ 
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SAY = 


= «i*( 


CAY + CAYBUR * PH. 


i) 


(33) 

(3U) 


T. = T. * exp (-SAY) 

1 i ^momo 


(35) 


where 


0 -^ 


= transmission for ith layer 

T, = transmission for all line-by-line calculations 

i,mono 

V = frequency, cm ^ 

CAYBUR = dummy variable 
SAY = dvimmy variable 

= ambient atmospheric pressure for layer i, atm 
2 

= molecules /cm of water vapor 
Az^ = vertical path length through ith layer 

s= cos 0^, where 0^ = angle through ith layer 

\ 

Wlj = water- vapor partial pressure, atm 


Program Description and Operating Instructions 

The CP program uses a Univac 1108 computer with four tape drives. The 
program tape is on unit C, the data tape is on unit A, and the scratch tape is 
on unit B. A tape drive Is also used to produce microfilm plot results on a 
General Dynamics 4060 plotter. If plot results are not desired or if the 
required hardware is not available, FORTRAN statements 335 to 348 should be 
removed. An option is available for producing resvilts on punched 
cards, if desired. 

The program is written in FORTRAN 5 language and, with storage, requires 
approximately 44K words. Only standard library routines compatible with 
FORTRAN 4 are used. No overlays are used. The subroutines MODATM, ATM0S3, 
INPUT, REFRAC, PATH, COSINV, SININV, Q, ALTITU, PRES, E,.R, and F are routines 
used to describe the atmospheric properties. These routines are described in 
reference Q.\ All these routines are usually used by MODATM, which is the sub- 
routine called by the main program CP. The purpose of each of these subpro- 
grams and its arguments are also defined in comment cards placed in the 
subroutine. 

The computer function XK (which is called by the main program CP) calcu- 
lates the Voigt spectrum line profile (eq. (2l)) by means of the routines XKl, 
XK2, and XK3. The arguments of XK, X and Y, are defined in equation (21). 

For detailed calculations involving the solar intensity reflection or 
the reflection of atmospheric emission as shown in equation (ll), multiple 
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and for upwelling intensity 





These functions are required for inputs to inversion routines that derive 
temperature structure from the upwelling radiance in the l4- to 15-micrometer 
wavelength region and water-vapor specific humidity structure from the upwell- 
ing radiance at wavelengths of 6.3 to 7.0 micrometers and a knowledge of the 
temperature structure. Chahine (ref. lO), Barnett (ref. ll), and Smith 
(ref. 12) have shown that it is necessary to have weighting functions for a 
class of atmospheric conditions, near the solution to ensure rapid convergence. 


Continuum 

The wave number region between 800 and 1250 per centimeter encompasses the 
continuum absorption due to water vapor. Calfee has used the continuum absorp- 
tion coefficient given by Burch (ref. 13) for this region. First, he uses the 

2 

optical mass of water vapor (U^ in molecules /cm ) to calculate the partial 
pressure due to water vapor (Wl^) for the ith layer. 



U^Tu. 


* 7.349 >< 10 
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(31) 


PH^ *= -Wl^ + 0.005 * + Wl^^ 


CAYBUR 


= -0.0665335 10“^^ + 0.3721108 


•X 10"^^v^ 


- 0.5971389 X 10 “^^ V + 0.3113485 10" 


20 


SAY = « (cay + CAYBUR PH 


i) 


^ * exp (-SAY) 

1 i,mono 


(32) 


(33) 

( 34 ) 

(35) 
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vhere = transmission for ith layer 

= transmission for all line-by-line calculations 

i,mono 

V = frequency, cm ^ 

CAYBUR = dummy variable 
SAY = dimimy variable 

= ambient atmospheric pressure for layer i, atm 
2 

= molecules/cm of water vapor 
AZj = vertical path length through ith layer 

= cos 6^, where 6^ = angle through ith layer 
Wl^ = water- vapor partial pressvire, atm 


Program Description and Operating Instructions 

The CP program uses a Univac 1108 computer with four tape drives. The 
program tape is on unit C, the data tape is on unit A, and the scratch tape is 
on unit B. A tape drive is also used to produce microfilm plot results on a 
General Dynamics it060 plotter. If plot results are not desired or if the 
required hardware is not available , FORTRAIJ statements 335 to 3U8 should be 
removed. An option is available for producing results on punched 
cards, if desired. 

The program is written in FORTRAM 5 language and, with storage, requires 
approximately U4K words. Only standard library routines compatible with 
FORTRAN are used. No overlays are used. The subroutines MODATIi* ATM0S3, 
INPUT, REFRAC, PATH, CQSINV, SINIFV, Q, ALTITU, PRES, E,.R, and F are routines 
used to describe the atmospheric properties. These routines are described in 
reference 8. All these routines are usually used, by MODATM, which is the sub- 
routine called by the main program CP. The purpose of each of these subpror- 
grams and its arguments are also defined in comment cards placed in the 
subroutine. 

The computer function XK (which is called by the main program CP) calcu- 
lates the Voigt spectrum line profile (eq. (2l)) by means of the routines XKl, 
XK2, and XK3. The arguments of XK, X and Y, are defined in equation (21). 

' For detailed calc\ilations involving the solar intensity reflection or 
the reflection of atmospheric emission as shown in equation (11), multiple 
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executions of the program are necessary. A separate run is required to calcu- 
late equations (l4) to (l6) for each required path through the atmosphere. 

The triangular instrument function is adequate for narrow-hand spectrom- 
eters; for instruments of wider hand pass, it is necessary to place a P in 
column 1 of card type 1. This step results in punching of two types of cards. 

The first type will punch v (cm”^), transmission, and transmitted radiance 

IX, F 9 . 2 , EI 3 . 3 , and EI 3 . 3 . The second type will punch v (cm ^), downwelling- 
’emitted radiance, and upwelling-emitted radiance IX, F9.2, E13.3, and E13.3. 
These cards can then he used with the appropriate filter functions for calcu- 
lations of the wide-hand response of instruments to various concentrations of 
gases and various model atmospheres. 

A source code listing of all required software is given in appendix A. 

The listing comprises 1355 FORTRAN cards. A detailed flow chart of the CP 
program is shown in figure 1. The number of inpiit cards required varies from 
as few as 20 that might he used in a 10- layer prediction of laboratory trans- 
mission for one spectral region ( 516 O to 5200 cm (appendix B, section 2) to 
approximately 35 that might he encountered when using a 10-layer prediction of 
atmospheric transmission through a radiosonde-specified atmosphere for' three 

Spectral regions (1150 to 1300, 550 to TOO, and TOO to 85 O cm (appendix B, 
section 3). Sample runs are provided in these same sections of appendix B for 
the input cards shown, and the generalized input format is presented in 
appendix B, section 1. 

A small section of the 15 250 lines of data stored on the data tape is 
presented in appendix C. A typical run of a 10-layer model using radiosonde 

data for three 150-cm~^ spectral regions will require approximately 10 minutes 
of computer time. However, run time is a function of the number of spectral 
lines in the chosen spectral interval, and some runs have required as many as 
U 5 minutes on the Univac IIO 8 computer. The number of pages usually produced 
is approximately 50 per spectral interval chosen, provided weighting functions 
are printed. If weighting functions are not printed out, approximately 10 
pages would he printed. Because the program is cyclic and returns to read 
additional spectral interval cards , the terminating conditions usually result 
from an attempt to read an end-of-file card. 

COMPARISONS TO TEST DATA AND LIMITATIONS 

OF THE MODEL ; ; 


Comparisons were made with laboratory spectra to perform necessary hut 
insufficient tests of the model. The model reproduces Calfee original water- 
transmission data (ref. 3) well for 0.01 and 0.1 centimeter of precipitahle 
water at 1 atmosphere pressure at frequencies between 5182 and 5193 per centi- 
meter (figs. 2 and 3). For low pressure, the model underestimated the ozone 
absorption by a factor of approximately 2 at wave numbers between 9^0 and 
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1120 per centimeter (fig. 1|). The CP model reproduces carton dioxide data veil 
for pressures of 1 . 00 , 0.085T» and 0.02053 atmosphere except at vave numbers 

between 600 and 660 cm”" for the lower pressures (figs. 5 to T)* One set of 

carbon dioxide lines at wave numbers between 6 U 0 and 650 cm” is evidently mis- 
sing from the data in figure 7* The comparison of spectra taken from an RB-5TF 
aircraft flying over the Gulf of Mexico to CP predictions from the Lake 
Charles, Louisiana, radiosonde on April 30, 1970, is shown in figure 8 . A com- 
parison of Skylab S191 spectra of Monroe Reservoir to those predicted by the CP 
program using radiosonde data from nearby Salem, Illinois, for June 10, 1973, 
is shown in figure 9* Because the amount of ozone in the atmosphere is unknown 
test cases were run for optical masses of both O.OlUU and O.lHU atm* cm. 

The results of Kunde and others (ref. 19) afford a comparison of a line- 
by-line model to Nimbus 4 interferometer data of 5 percent in the U25“ to 550- 
per-centimeter water continuum and the 750 - to 1200 -per-ceivcimeter atmospheric 
window and 5 to 10 percent in the 667 -per-centimeter carbon dioxide absorption 
region. The absolute accuracy was 5 to 10 percent for each of the parameters: 
measured radiances, in situ measurements, and atmospheric transmittances . The 
conclusion was that it is not possible to uniquely determine the error arising 
from each source and that it is, therefore, very difficult to improve atmo- 
spheric transmittances (predicted by a computer program) through the radiance 
comparison technique. The results shown in figure 9 indicate a comparison 
within 5 percent for the atmospheric wixidow, but reveal larger errors outside 
this region. Whether these errors arise from the CP program' or the S191 spec- 
trometer or from both Is unknown. Some off-band contamination of the 6 - to 9“ 
and 13.0- to 15 . 5 -micrometer regions of the SI 9 I data has been identified that 
causes radiance of these regions to be excessively high. The amount of this 
correction is currently being determined by the sensor performance personnel at 
the NASA lyndon B. Johnson Space Center. 


CONCLUDING REMARKS 


‘ The results of testing the CP program indicate that the best results can 
be obtained in the near-infrared water bands. The absorption due to -water va- 
por and carbon dioxide in the thermal infrared band appeared to be less reli- 
able in comparison to spacecraft-acquired data and band models; however, nei- 
ther of these tests is sufficient. Comparisons of laboratory carbon dioxide 
transmission in the thermal infrared band show good agreement except in regions 
where lines are known to be missing. The comparison of ozone transmission at a 
wavelength of 9*6 micrometers to laboratory data showed unexpectedly large 
differences . 

No testing has been done for the nitrous oxide or methane transmissions. 

At present, the data used in the program are probably better known for carbon 
dioxide, nitrous oxide, and water vapor than for ozone and methane. Although 
tests for all molecules and spectral regions have not been performed, the tests 
have been sufficient to ensure that the program mechanics_ are sound and work 
well on the Univac IIO 8 computer. Conversion to FORTRAN U for CDC or IBM 


3 ^ 


computers will no doubt require new tape read-and-write software. However, the 
necessary changes should be simple, few, and obvious, because most FORTRAN coding 
usied here did not involve sophisticated FORTRAN 5 logic. 

Lyndon B. Johnson Space Center 

National Aeronautics and Space Administration 
Houston, Texas, December 27, 197^ 

951-16-00-00-72 
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TABLE I.- SUMMARY OF UNIT RELATIONSHIPS 


1 precipitable centimeter water = 1 g*cn w 


Symbol and previous unit 


Multiplying factor 
(b) 


Current unit 


U, atm* cm STP 


1.219 X 10"^ X ~ 
o 


U, g/cm*^ 


U, atm* cm STP 


82.06 X ^ 


atm* cm STP 


2.689 X 10^^ molecule* pm 


2 

U, g/cm water 


3.3*» X 10^^ 


molecule* cm 


Kq(v), (atm* cm)” STP 


82.06 5 - 

M 


(g*cm"^) 


Kq(v), (g*cm”^) 


, 1.219 X 10 


-2 M 


(atm*cm)”^ STP 


Kq(v), (atm^cm)”"^ STP 


3.72 X 10 


^2 

(molecule* cm ) 


Kq(v), (g*cm“^) 


(molecule* cm ) 


Kq(v), (atm*cm)"^ STP 356.3 | 


dB/(g*cm” ) 


Kq(v), (atm*cm)” STP 


U. 3 U 3 


dB/(atm*cm) STP 


Kq(v), (g*cm“^) 


U. 3 U 3 


^p 

dB/(g*cm” ) 


S, STP 

atm* cm 


82.06 5 - 


1.219 X 10 


-2 M 


atm*em 


S, STP 

atm* cm 


3.72 X 10 


molecule* cm 


molecule* cm 


uata from reference 5. 

b 

9 = temperature, K, M = molecular weight, A * Avogadro's number 
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SlNAL PAGE IS POOR 









Read in spectral regions 
and specifications 


Read in target temperature, 
emissivity location , and 
sensor location 


Calculate path length in 
each atmospheric layer 


Target 

and sensor altitudes 

s. <0.00001 km ^ 


Read in gas concentrations 
for each layer and calculate 
gas mass in each layer 


Calculate continuum 
800 to 1250 cm ^ 


Read in spectral region 
and specifications 


Calculate Lorentz or Voigt 
broadening and calculate 
absorption coefficient for 
each frequency 


Calculate monochromatic 
transmission through 
entire atmosphere 


Calculate atmospheric 
emission terms 


Integrate over triangular 
filter function 


Print out, punch, and 
plot (microfilm) results 


Figure 1.- Flow chart of the CP program. 











O CP program calculation 



Frequency, cm 


Figure 2.- Comparison of the predicted absorption for program CP with the original Deutschmann 
and Calfee (ref. 3) single-layer water-vapor-ahsorption spectrum for the frequency range 
5182 to 5193 per centimeter at a pressure of 1.0 atmosphere, at a concentration of 0.01 
centimeter of precipitable water vapor for a temperatxire of 287. T K, and for a spectral 


slit width 


of 0.5 cm' 




O CP program calculation 



0 *— 
5182 


5184 


Frequency, cm 


Figure 3.- Comparison of the predicted ahsorption for program CP vith the original Deutschmann 
and Calfee (ref . 3) single-layer water -vapor-ahsorption spectrum for the frequency range 
5182 to 5193 per centimeter at a pressTore of 1,0 atmosphere, at a concentration of 0.1 
centimeter of precipitable water vapor for a temperature of 287.7 K, and for a spectral 

slit width a of 0,5 cm 





Transmittance 


Experimental 
Variable a 
a = 0.08 cm ^ 
a = 0.06 cm ^ 

CP program calculation 


Wk 


Wave number, cm 


Figure 5.- Comparison of predicted transmission for program CP vith the carton dioxide 

transmission of Drayson and others (ref. l6) for an eq^ui valent pressure of 1.00 atmosphere 
and an optical mass of 106.2 atm»cm_,._ „. 





Transmittance 


— Experimental 
— — Variable a 
■ a = 0.08 cm ^ 

a = 0.06 cm ^ 

O CP program calculation 




Wave number , cm 

Figure 6.- Comparison of predicted transmission for program CP with the carbon dioxide 

transmission of Drayson and others (ref. l6) for an equivalent pressure of 0. 08 57 > atmosphere 
and an optical mass of 212.1 atra*cm_„^ 
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dioxide band missing 
from CP table of lines 
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Wave number, cm 


— Experimental 

— Variable a 

■“ «= 0.08 cm ^ 

••• a = 0.06 cm ^ 

CP program calculation 


Figure 7«- Comparison of low-resolution spectra (ref. 17) with theoretical calcxilations of 
Drayson and others (ref. l6) and CP program for an equivalent pressure of 0.02053 atmosphere 
and an optical mass of 6.30 atm-cm_-^ 
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O Intensity from the filter-wheel spectrometer obtained by using preflight 
and postflight calibrations 


— — Nonlinear least-squares-fit black-body temperature - 297.1966 K 

- — Predicted intensity upwelling from the Gulf of Mexico (T = 298.5 K) 

accounting for transmission and emission of the atmosphere using radiosonde 
data from Lake Charles/ Louisiana, April 30, 1970, at 00:00 G.m.t. 

Figure 8.- Cali'bration performed on April 30, 19T0, using the Gulf of Mexico 
as a sotirce and accounting for the intervening atmospheric transmission 
and emission. 
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Figure 9.- Comparison of S191 sensor output data for Monroe Reservoir vith CP model and 
Anding model radiosonde data from Salem, Illinois, on June 10, 1973. 
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MEASUREMENT OF OCEAN PARAMETERS IN SUPPORT OF SKYLAB MISSIONS 
IN THE GULF OF MEXICO OFF THE TEXAS COAST 

C. Laurence Kerb,* and John F, Potter** 

INTRODUCTION 

A series of four separate ocean truth missions were 
conducted to obtain the ground truth data necessary to 
support this investigation. The first mission was conducted 
on 5/31/73 from Freeport, Texas along a line extending 
southeasterly from Freeport, Texas at 28®55'N, 95®18'W 
to approximately 15 miles offshore at a 28®45'N, 95°09'W. 

The data obtained included measurements of chlorophyll 
concentration, turbidity, water temperature, water samples, 
and an attempt to measure incident spectral solar irradiance. 

Skylab data was not obtained for this mission and no laboratory 
’spectral analysis was performed on the water samples which 
were taken. 

The second mission was conducted on 6/5/73 from Port 
O’Connor, Texas and consisted of data taken in Matagorda Bay 
.in the vicinity of 28®27'30"N, 96°21’30''W and coastal data 
taken along a line extending nearly due east to approximately 
2^®23*40"N, 95®55'30”W. The data taken included measurements 
of chlorophyll concentration, turbidity, water temperature, 
water samples, equivalent Lambert reflectance of the ocean, 

and solar photometer measurements. 

* Goddard Space Flight Center 
**Lockheed Electronics Corporation 
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A Skylab data -take was obtained for this mission. 
Unfortunately, the data take was at a solar elevation angle 
of 82® and thus the data was badly contaminated by solar 
glitter from the ocean. Extreme problems were also 
encountered in attempts to measure the spectral reflectance 
of the ocean from shipboard due to specular reflection 

by the ocean. The difficulties of these reflectance 

measurements were further compounded by ship motion. 

The third mission was conducted on 8/8/73 from Freeport, 

Texas along the same line that was used for the first 
mission. The data obtained was limited to collecting water 
samples and measurements of water temperature since the vessel 
being used for data collection experienced a total power 
failure early in the mission, rendering all auxilliary 
instrumentation that required power unuseable. After 
returning to dockside later in the day, however, measurements 
of chlorophyll concentration and turbidity were performed 
using the water samples that had been collected. The 
weather during this third mission was extremely poor. The 
weather was dominated by thunderstorms early in the morning, 
and scattered rains and heavy cloud cover during the time 
ofj the Skylab overpass. A Skylab data -take was obtained for 
this mission, but the data was of no value for oceanography 
due to the presence of the heavy cloud cover. 

The fourth ocean truth mission was conducted on 9/16/73 
in the vicinity of 120 nm off the coast of Port O'Connor 
Texas along a line nearly due west from 27®15'N, 94®35'W 
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to 27®19'N, 94®54'W, and then along a second line extending 
northeast to 28°04'N, 95®52'W. The data obtained from this 
mission included measurements of chlorophyll concentration, 
turbidity, water temperature, water samples, equivalent 
Lambert reflectance of the ocean, solar photometer measurements, 
arid extinction of the ocean over a one meter path at various 
depths. The v/eather for this mission was excellent with 
clear skies, winds in the 5 mph range, and a solar 
elevation angle of 45® at the time of the scheduled Sky lab 
overpass at 10:45. Unfortunately, a Skylab data-take was 
not performed due to a weather forecast which predicted poor 
weather for the test site area. 

• i •' 

. i 
i 

: t: 

RESULTS AND DISCUSSION 

I The data acquired on Missions 2 and 4 covers the range of 

; ‘ "ocean parameters encountered in bay, coastal, and offshore 

I sbudies and is presented in this section. 

Table I gives the latitude, longitude, and time 
associated with each site for. which measurements were made 
in the vicinity of Matagorda Bay and the associated coastal 
offshore waters on June 5, 1973. 

Table II gives the measured values of turbidity, chlorphyll 
concentration, and water temperature which were found at these 
sites. The value of chlorophyll concentration were measured 



TABLE 1 

SITE NUMBER, LONGITUDE, LATITUDE, AND TIME FOR GROUND TRUTH DATA 
COLLECTION ON JUNE 5, 1973, IN THE VICINITY OF MATAGQRDA BAY 


SITE 

LONGITUDE 

LATITUDE 

TIME 


96°21’40" 

28°27’40" 

9.15 

2 

96°21'30" 

28°27'40" 

9.20 

3 

96°21’30" 

28°27'30" 

9.35 


96°20'40" 

28°26’20" 

10.15 

5 

96°19'30" 

28°25'00" 

10.25 

6 

96°19’00" 

28°24'20" 

10.30 

7 

96°15'00" 

28°24'20" 

10.50 

8 

96°l3'00" 

28°24'20" 

11.10 

9 

96°11’10" 

28°24'20" 

11.30 

10 

96° 8 '10" 

28°24'10" 

11.50 

11 

96° I’lO" 

28°23'50" 

12.20 

12 

95°55'30" 

28°23'40" 

14.20 

13 

95°56’50" 

28°24'20" 

14.30 

14 

95°58'50" 

28°25'20" 

14.44 

15 

96°00'30" 

28°26'20" 

14.55 

16 

96°02'10" 

28°27’10" 

15^3 

17 

96°04'10" 

28°28'00" 

16.00 

18 

96°09'30" 

28;°28'40" 

16125 

19 

96°11'50" 

28°27'20" 

16 1 30 

20 

96°13'00" 

2$°27'40" 

16.40 

21 

96°14’30" 

28°26'50" 

16.55 

22 

96°15'10" 

28°26’40" 

17.05 

23 

96°15’10" 

28°26’40" 

17,10 

24 

96°17'10" 

28°25’30" 

17.25 

25 

96°17'50" 

28°23’00" 

17.35 

26 

96®19'30" 

28°24'50" 

17;. 50 

27 

96°19'30" 

28°25'00" 

17.55 

28 

96°20’10" 

28°25’50" 

IQ.OO 

29 

96°21'40" 

28°27'00" 

18.10 

30 

96°20'50" 

28°27'50" 

18.20 
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TABLE II 


TURBIDITY, CHLOROPHYLL, AND TEMPERATURE OF OCEAN GROUND TRUTH 



SITES IN THE VICINITY 

OF MATAGORDA BAY, JUNE 

5, 1973 


SITE 

TURBIDITY ( JTU) 

CHLOROPHYLL 


TEMPERATURE 

1 

17.0 

6.9 



2 

19.0 

5.0 


- 

3 

13.0 

6.9 


- 

4 

19.0 

7.5 


27.7 

5 

23.0 

6.2 


- 

6 

16.0 

7.1 


27.7 

7 

12.0 

5.9 


27.5 . 

8 

5.5 

1^2 


- 

;9 

5.0 

- 


26^9 

lo 

4.0 

2.5 


26l7 

11 

2.4 

1.0 


26.7 

12 

0.9 

1;2 


27.1 

13 

1.4 

1.8 


- 

14 

1.7 

1;8 


- 

15 

1.5 

2.4 


26.9 

16 

3.5 

4.0 


26.9 

17^ 

2.6 

3.9 


27.1 

18 

5.5 

5.4 


27.7 

19 

6.5 

5.4 


29. 2 

20 

25 

7.7 


28.3 

21 

22 

13.0 


28.3 

22 

20 

14.5 


28.4 ' 

23 

7.5 

9.0 


28.4 

24! 

6.5 

9.0 


28.4 

25 

I ■ 6.5 

6.0 



26 

I6i' 

8.2 


28.4 

27 

15: 

8.2 


28.3 

; j ■ ■ 

28i 

11 

7.1 


28.3 

29 

9.5 

6.9 


28.3 

30 

9.0 

11.2 


28.6 


using a Turner fluorometer. The fluorometer was calibrated 
periodically during the course of each mission using distilled 
water to determine zero offset readings. The relation between 
chlorophyll concentration and calibrated fluorometer readings 
was established by laboratory analysis of water samples by 
Texas A&M University. The measurements of turbidity were 
likewise made in situ using a Hach Laboratory Turbidemeter , 
Model 2100, which measured the intensity of light scattered 
by a sample of ocean water at an angle of 90° relative to 
the direction of evidence. The turbidity readings were 
calibrated relative to a sample of Formazin solution and 
the calibrated measurements are given in Jackson Turbidity 
Units (JTU) . The measurements of water temperature were 
made using a calibrated thermometer and the bucket method. 

As may be seen from Table II, a values of chlorophyll 
found were in the range of 1.0 to 14.5 mg/m which indicates 
a very high level of phytoplankton activity. The associated 
values of turbidity range from as low as 0.9 to as high as 

2^.0 JTU. The range of ocean temperatures found varied from a 

1 

low of 26.9° to a high of 29.2°C. It may be noted, however, 
that the effect of solar heating during the day accounted 
for about a 0.6°C change in ocean temperature. 

It may also be seen from Table II that trends occur 
in both the levels of chlorophyll and turbidity. The levels 
of chlorophyll and turbidity associated with sites 1 thru 7 
are very high. These sites start in Matagorda Bay and move 


7 


5 ^ 

into the waters adjacent to Matagorda Peninsula. The values 
of chlorophyll and turbidity then decrease as the sites move 
further offshore (sites 8 thru 11) and a minimum value of 
turbidity is found at site 12 which is the farthest point 
from landfall. The values of chlorophyll and turbidity then 
slowly increase as one approaches site 17 (4.5nm from the 
Peninsula ) , and then have a large increase in value in the 
vicinity of site 20 through 22 which lie only about 1.5nm 
off the coast of the Peninsula. 

The high levels of turbidity which are associated with 
Matagorda Bay and the coastal waters within 2nm of the coast 

! of the Peninsula may be understood as naturally occuring in 

I . 

’ shallow waters due to runoff from the shore and also 

: perhaps due to mud and silt. It is also important to note, 

► 

, however, that high values of turbidity and high values of 

chlorophyll have a tendency to go together. The relationship 
between turbidity and chlorophyll was investigated by plotting 
(Figure 1) the values of turbidity vs. chlorophyll for each 

« o|f the data sites in Table II. When a correlation coefficient 

• . --r ... ..... .. 

was calculated for this data, it was found to be 0.68. This 
value indicates a correlation between the data significant at 
a confidence level of greater than 99%. 

One would not a priori hypothesize a relationship such 
as that found in Table II and shown in Figure 1. It is 
reasonable on a physical basis that such a relationship 
might occur, however, since high or low turbidity is 



1 



FIGURE 1 - Chlorophyll vs. turbidity 
for Matagorda Bay and . 
associated coastal water. 


„ Turbidity (JTU) 


associated with correspondingly high or low values of 
nutrients in the coastal watets. The high level of 
nutrients provide favorable conditions for high rates of 
photosynthesis and thus of chlorophyll. 

Water samples were collected at each site during mission 
2 and these samples were returned for laboratory spectral 
analysis. This analysis consisted of measuring the spectral 
transmittance of a 10 cm sample of ocean water using a 
Cary-14 spectroradiometer . The results of these laboratory 
measurements are shown in Figures Al through A27, In , 
order to eliminate the effects of absorption by the 
cell windows, reflection losses at the cell interface and 
absorption due to the molecular absorption bands of liquid 
water, the transmission measurements were made relative to 
a reference cell containing a 10 cm path of distilled water. 

As is shown in Figures Al through A27 , the resultant 
spectra generally show a nearly linear relationship between 
transmission and wavelength over the spectral region from 
4000 to 8000 A® with higher values of transmittance at the 
longer wavelengths. The prime feature of these spectral 
curves is the shift in the overall level of transmitted 
energy between sites. At a wavelength of 4500A® the 
transmission varies from a low value of about 18% for site 2 
in Matagorda Bay to a high value of 95% for the offshore 
coastal waters. The overall shape of the transmitted 
spectrum is consistent with the effects of large particles 
Mie scattering. 


Figures A1 through A27 

The ratio of the energy transmitted by a 10 cm sample of 
ocean water to that transmitted by an equal path length of 
distilled water as a function of wavelength for the spectral 
region from 4000 to 8000 A®. 
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The results of mission 4 are given in Tables III and IV 
and in Figures B1 through B9. This mission was conducted 
on 9/16/73 approximately 120nm off the coast of Port O'Connor, 
Texas, Table III gives the position coordinates and time 
associated with each site. Table IV gives the measured 
Values of chlorophyll concentration, turbidity, and 
water temperature which were found at each site. 

As is shown in Table IV, the values of chlorophyll 
concentration were in the range of 0.16 to 0.48 mg/m^ for 
the more distant sites 1 through 4 , and a value as high as 
1.53 mg/m^ was found at site 5 which was measured about 
35 nm out of Port O'Connor on the return trip. The values 

' I ' 

of turbidity found were in the range of 0.15 to 0.50 JTU. 

The range of ocean temperatures measured varied from 28.6°C 
to 29.4®C. 

A comparison of the chlorophyll and turbidity levels found 
in the offshore waters of mission 4 (Table IV) with those 
found in the coastal and bay waters of mission 2 (Table II) 
indicates that the levels of chlorophyll and turbidity 
occuring in the offshore waters of the Gulf of Mexico are 
between one and two orders of magnitude lower than those - 
found in the coastal waters. Also, it may be observed 
from Table IV, that there is no evident correlation of 
turbidity and chlorophyll values for the offshore regions 
as was found in the coastal regions. 
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TABLE III 

SITE NUMBER, LATITUDE, LONGITUDE AND TIME FOR GROUND TRUTH 
DATA COLLECTION ON SEPTEMBER 16, 1973, IN THE VICINITY OF 
120 nm OFF THE COAST OF PORT O’CONNOR, TX 


SITE 

LATITUDE 

LONGITUDE 

TIME 

1 

27°15' 

94°35' 

8.10 

2 

27°10' 

94°41’ 

10.10 

3 

27°09' 

94:°42 ’ 

10.42 

4 

27®19' 

94° 54' 

12.20 

5 

28°04' 

95°52' 

16.50 


% 



I 


I 


J 
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TABLE IV 


TURBIDITY, CHLOROPHYLL, AND TEMPERATURE OF OCEAN GROUND TRUTH 
SITES ON SEPTEMBER 16, 1973, IN THE VICINITY OF 



120 ntn OFF TR 

M COAST 

OF PORT O'CONNOR, TX 


SITE 

TURBIDITY 

(JTU) 

CHLOROPHYLL (ms/m^) 

TEMPERATURE 

1 

0.47 


0.16 

28.6 

2 

0.50 


0.36 


3 

0.15 


0.48 

29.3 

4 

0.32 


0.24 

29.4 


5 


0.17 


1.53 


29.3 



,1 .1 I 1 
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Spectral transmission curves of the water samples which 
were collected during mission 4 are shown in Figures B1 
through B9. For these laboratory measurements, the total 
transmission of the water samples were measured rather than 
making a ratio measurement as was done previously. The 
spectrum of distilled water is shown in Figure B1 for 
comparison purposes. As is shown in these figures, the 
molecular absorption spectrum of water dominates the overall 
spectrum with relatively high transmission in the visible 
portion of the spectrum, 4000 to 6500A°, and strong 
absorption features in the near infrared located near 
7600A®, 9800A°, and beyond 11,500A°. For most sites, 
two water samples were taken and the spectral curves for 
different samples from the same site generally show agree- 
ment to within 5%. A spectrum of a water sample which 
contained kelp is shown in Figure B6 . . As may be seen from 
this figure, the spectrum is dominated by strong absorption 
_in the blue and green portions of the spectrum. 

A problem encountered in making laboratory analyses 
of Water samples which contain living organisms is the 
time delay between acquiring the water samples and per- 
forming the laboratory analyses. In the case of the 
spectral analsis which we performed, this time delay 
was between one and two weeks. In order to preserve the 
water sa; pies, they were refrigerated during the return 
trip and then frozen until the laboratory measurements 
could be started. It is likely, however, that changes 



Figures B1 through B9 


The transmission of a 10 cm sample of ocean water as a 
function of wavelength for the spectral region from 4000 to 
10,000 A® including effects of cell transmittance. 




BliPKODUCIBILITY OF THE 
OBIGWAL PAGE IS POOR 


0 5 


0.4 


0 3 


02 


0.1 


on 


pic 32 

fcATE lO-l-TB TIME 


PROJE CT tOIOBg, Sa/A^E BlACP- 
RANGE . 4 - 1 . 2 a , AN'GL E 


SCAN RATE PCRie» ^ 

SLIT CONTROL -ZS OPclUT»R N -C . 

■PREPARED BY: 

MEASUREMENT ACQUISITION LABOR^VIORY 
DATA APPUCATIONS f, PHYSICS BRilNCH, EOD 
NASA, LYNDON B. JOHNSON SPACE CENTER 
HOUSTON, TEXAS 







Ft(. 9 y' ■ ■ 

0 .4 5 !/» ~ I ^ ^ _ 0 4 

; DAt^ l Q-l-73 TIM: _ ~ ~~~ 

rnoJECT j^A3E&_samplej^^^ST ^ i 

Range -4-I.Z A: angle —::: •• 

^ , SCAN Bate A/s ec- PEK io>_r I 

— SLIT r.CNTROL 2-^ - Or'L.BArek „_J 



PREPARiiD by: 

0 2 MEASURi-Ml-NT ACQUISITION LABORATORY 
- 7-- UATA APPLICATIONS 5 PHYSICS BRjVNCII, BOD 
NASA, LYNDON B. JOHNSON SPACE CENTER 
HOUSTON, TEXAS 


[ 








Cn.iC.v».- Tn. , . . CHART No Hco 

li 





0 



nrpRonuaBiuTY of thb 

OWOINM- PAGE IS POOR 


4 


?/ 




ivj 


i 


1 


1 


i 

fv/ 


... « ... . CH/<M »v I4gf) 




0 5 0 5 

' ficey ■ ' ' 

DATt 10>|-~?3 TIME 

P.A tf/-T_tijJkTE^e. y.AM-’tt Ct-gOMT 0 4 

R ANGELA — 

SCa/^ »Are 25^EC.»’l*i-)» .? 

sm control anRATjR K) c.-_^ 

03 03 

PREPARED BY: 1 *' ~ 

MEASUREMENT ACQUISITION LABORATORY 

DATA APPLICATIONS 6 PHYSICS BRANCH, EOD ‘ 

NASA, LYNDON B. JOHNSON SPACE CENTER 
0 2 HOUSTON, TEXAS 0 2 





ORIGINAL PAGE IS 
OF POOR QUALITY 







CMAm K. 14,^ 

tl 



I 


0 *; 

Ft<* (fV 

MATf^i tY -'»//<- 

DATt \0-\~m TI*A! 0 4 

RANGt . 4-1 - Z-A *NGit 

SCAN K*U 

SUt CONT*#V_'2iS__OPt8AlOR..W;£_l, 
03 


PREPAREP BY: 

MEASUREMENT ACQUISITION LABOR.\TORY 
DATA APPLICATIONS f, PHYSICS BRANCH, EOD 

- NASA, LYNDON B. JOHNSON SPACE CENTER 

HOUSTON, TEXAS 


0 t 


; t 




r\ 

- 1 ■ \ 
/ 



0 •> 


c: 


o I 








I J 



-52- 

tock place in the biological composition of the water samples 
during the holding period and thus the measured spectral 
transmitances are not the same as would have been observed 
if the measurements had been made in situ. 

A considerable effort was made during the^ course of 
this investigation to develop a technique for measuring 
the scattering properties of the water samples. In order 
to perform these measurements , a Science Spectriam Light 
Scattering Photometer was rented from the manufacturer. 

The instrument was capable of making qualitative measure- 
ments of the shape of the scattering properties of liquid 
samples over the angular range from approximately 1 to 170°. 
A number of difficulties were encountered, however, in 
attempts to obtain quantitative data. These problems 
were in the areas of instrument noise which generated 
lafge spurious peaks in the scattering measurements; 
highly, nonlinear amplification electronics which required 
the development of a new calibration technique in order 
to compensate for this effect; and, spurious scattering 
•which was caused by the introduction of neutral density 
filters which were needed to attenuate the intense signal 
which occurs for small scattering angles. The primary 
instrument problems were fixed and calibration techniques 
for using the light scattering photometer were developed 
over the period from August to December, 1973. However, by 
the time the photometer could be put into operation, the 



53 




water samples which had been acquired over the course of the 
missions had deteriorated. 


CONCLUSIONS 

High values of chlorophyll and turbidity were measured in 
the vicinity of Matagorda Bay# Texas, and the associated 
offshore coastal waters. Values of chlorophyll were found 
to be in the range of 1.0 to 14.5 mg/m and the values of 
turbidity were in the range from 0.9 to 25.0 JTU. A 
correlation coefficient was calculated for the chlorophyll 
and turbidity data which yielded a value of 0.68. This 
indicates that the levels of chlorophyll and turbidity 
were correlated at a confidence level of greater than 99% 
for the bay and near coastal data. This result, if it is 
generally characteristic of coastal regions , would have 
important implications for studies concerned with the remote 
sensing of chlorophyll in these regions. It may be noted 
that high levels of turbidity can be identified easily from 
remote sensing data. Therefore, if a relationship between 
turbidity and chlorophyll exists such as was found in our data, 
then the presence of high turbidity levels in remote sensing 
data could be used to infer the presence of correspondingly 
high values of chlorophyll. 

Measured values of chlorophyll and turbidity were also 
obtained for the offshore region approximately 120 nm off 



the coast of Port O'Connor. For this distant offshore 
region, the values of chlorophyll were in range from 

3 

0.16 to 0.48 mg/m and the values of turbidity were in 
the range from 0.15 to 0.50 JTU. For the relatively 
low values of chlorophyll and turbidity occuring in these 
offshore waters, there was no significant correlation 
between the values of chlorophyll and turbidity. 

For the offshore area and also for the coastal data 
previously discussed, sea surface temperatures were 
measured and found to lie within a) narrow ranges of less 
than 2^C. No important correlation was fo;ond between 
sea surface temperature and chlorophyll . 



